On the Hydrodynamic Interplay Between a Young Nuclear 
Starburst and a Central Super Massive Black Hole 



Filiberto Hueyotl-Zahuantitla 1 ,Guillermo Tenorio-Tagle 1 ' 2 , Richard Wunsch 3 , Sergiy Silich 1 

& Jan Palous 3 



ABSTRACT 

We present ID numerical simulations, which consider the effects of radiative 
cooling and gravity on the hydrodynamics of the matter reinserted by stellar 
winds and supernovae within young nuclear starbursts with a central supermas- 
sive black hole (SMBH). The simulations confirm our previous semi-analytic 
results for low energetic starbursts, evolving in a quasi-adiabatic regime, and 
extend them to more powerful starbursts evolving in the catastrophic cooling 
regime. The simulations show a bimodal hydrodynamic solution in all cases. 
They present a quasi-stationary accretion flow onto the black hole, defined by the 
matter reinserted by massive stars within the stagnation volume and a stationary 
starburst wind, driven by the high thermal pressure acquired in the region be- 
tween the stagnation and the starburst radii. In the catastrophic cooling regime, 
the stagnation radius rapidly approaches the surface of the starburst region, as 
one considers more massive starbursts. This leads to larger accretion rates onto 
the SMBH and concurrently to powerful winds able to inhibit interstellar matter 
from approaching the nuclear starburst. Our self-consistent model thus estab- 
lishes a direct physical link between the SMBH accretion rate and the nuclear 
star formation activity of the host galaxy and provides a good upper limit to the 
accretion rate onto the central black hole. 

Subject headings: galaxies: active — galaxies: starbursts — accretion — hydro- 
dynamics 
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Introduction 



Powerful starbursts have been conclusively detec ted in the nuc l ear re gions of galax- 



i es wi th active galactic nuclei (AGN; see for a revie w IVeilleux et al.l (120051) and [H eckman 



2007 



Watabe et al 



formation rate (Alexander et al. 



2008; 



( iCid Fernandes fc Terlevich 



1995 



2005 



(20091 )') . These include qu a sars (IHao et al.ll2005l . 120081 ). Seyferts (Ilmanishill2003l ; iDavies et al. 



Chen et al. 2009), submil 



Walter et al. 



imeter galaxies with an extreme star 



20090 and even low luminosity AGNs 



Gonzalez Delgado et al.l 120041 ) . Some s upermassive blac k 



holes (SMBHs) coincide with massive and compact nuclear star clusters ( jSeth et al.l 120081 ). 
some of them with a co mplex history of star f ormation are found in a number of near by 
spiral ( jRossa et al.ll2006l ) and elliptical galaxies (jWehner fc Harris! l2006t ICote et al.l 120061 ) . 



The interplay between young nuclear starbursts and their central SMBHs is not well 
understood and remains one of the central issu e in the theory of AGN galaxies, and in 



the p rocess of cosmological growth of SMBHs (ILipari fc Terlevich 



4' 

2009) and their co-evo l ution with the bulges of their host galaxies (IBegelman fc Nath 



2006; 



Booth &; Schave 



2005 



Di Matteo et al.ll2005l ; ISomerville et al.ll2008l ). In fact, even basic issues regarding, for ex- 
ample, the impact of type II supernovae (SN) on the matter left over from star formation, 
seem to be still undecided. One can find in the literature massive starbursts with a SN 



rate of 1 yr _1 structuring a gaseous disk of just 10 7 M Q (IWada fc Normanll2002l ) while other 
calculations a ssume that type II supe rnova might evacuate most of the nuclear region from 
gas and dust (jSchartmann et al.l 120091 ) and have considered only the evolution after the end 
of the type II SN epoch. 



Silich et al.l ( 120081 ) presented a self-consistent, spherically symmetric stationary solution 



for the gaseous flow around a SMBH at the center of a young (< 40 Myr) starburst region. 
They also found a threshold mechanical luminosity (L cr i t ) which separates systems evolving 
in the quasi-adiabatic regime with a small stagnation radius, R s t (at which the velocity is 
zero km s -1 ) and consequently a low accretion rate, from those with a large R st and a high 
accretion rate whose hydrodynamics are strongly affected b y radiative cooling (see their Fig- 
ure 6). The self-consistent semi-analytic solution found by ISilich et al.l (120081 ) is only valid 
for nuclear starbursts with a mechanical luminosity Lnsb < L cr i t . Above the threshold lumi- 
nosity, strong radiative cooling drives the gas thermally unstable. This has a strong impact 
on the thermal pressure gradient and thus on the location of the stagnation radius, facts 
which inhibit a complete semi-analytic solution. Thus cases with a mechanical luminosity 
Lnsb > L crit demand of a numerical in t egrat ion of the flow equations. Nevertheless, from 
the semi-analytic results of ISilich et al.l (120081 ) it was suggested that the thermalization of 
the kinetic energy released within the volume occupied by a young nuclear starburst (NSB) 
should lead to a central well regulated accretion flow onto the massive object, while a power- 
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ful wind, capable of preventing the accretion of the ambient interstellar gas (ISM) onto the 
SMBH, should emanate from the outer regions. Here we follow these suggestions and mea- 
sure the power of the resultant winds, as well as the accretion rates onto the SMBH and show 



(2000 


), 


Ciotti et al. 


(2009) 



gas is not too large and here we present results from our ID numerical simulations as a 
first approximation. The simulations confirm Silich's et al. (2008) semi-analytic results and 
extend them to starbursts with a mechanical energy input rate, L NS b > L crit . In this paper 
the input physics account for the mechanical feedback that a young nuclear starburst may 
provide to the accretion flow while allowing for radiative cooling and the gravitational pull 
from both t he central SMBH and the nuclear starburst. The feedback p rovided by the cen- 



tral SMBH ( ISilk fc Reeslll998l : ICiotti fc Ostrikerll200ll ; ICiotti et al.ll2009h will be the subject 



of a forthcoming communication. 

The paper is organized as follows: the numerical scheme used and the initial and bound- 
ary conditions are presented in section 2. Section 3 deals with the results from the simulations 
and compares them with the semi-analytic results. In section 4 the accretion rates and lu- 
minosities of the central SMBHs are obtained, as well as a measure of the possible impact 
of the resultant winds onto the host galaxy ISM. Section 5 gives our conclusions. 



2. Numerical Approach 



The numerical model is based on the finite difference Eulerian hy dro dynamic code 
ZEUS3D v.3.4.2, which solves the set of the hydrodynamic equations (IStone fc Norman 
19921 : Iwiinsch et al.ll2008h : 



dp/dt + V • (pit) = q m , 

du/dt + (it ■ V)u + VP/p = -V<p BH+NSB , 

de/dt + V ■ (eu) + PVu = q e -Q, 



(1) 
(2) 
(3) 



where, q m and q e are the mass and energy deposition rates per unit volume, e is the internal 
energy, Q = rtjn e A(T, Z) is the cooling rate, rtj and n e are the ion and electron number densi- 
ties, and A(T, Z) is the cooling function, which depends on the thermalized gas tempera ture, 
T, and metallicity, Z. The Raymond & Cox cooling function tabulated by iPlewal (119951 ) has 
been used in all calculations. The right hand side of equation ([2]) represents the gravita- 
tional acceleration a g = —GM r /r 2 , where M r = M BH + M NSB (r / R NSB ) 3 is the mass within 
a volume of radius r. Outside of the starburst volume, r > Rnsb, it is M r = M BH + M n $b- 
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All simulations have been carried out in spherical coordinates with symmetry along the 
6 and (f) directions and with a uniform grid in the radial direction. The c alcula tions account 
for fast radiative cooling (see iTenorio-Tagle et al.l ( 120071 ); IWiinsch et al.l (120081 )) and include 
the gravitational pull from the central SMBH and from the starburst, assuming that the 
stars are homogeneously distributed within a spherically-symmetric volume. 



2.1. Initial and Boundary Conditions 



The initial distributions of velocity, pressure, temperature and d ensity for all calculations 
were taken from the semi- analytic wind solution (ISilich et al.l 120041 ). without accounting for 
the gravitational pull from the starburst and the central SMBH. The initial condition is 
adapted to starbursts of the required size (see Table 1), with a mechanical luminosity Lnsb < 
L crit . The mechanical luminosity Lnsb = (where, Va,oo is the adiabatic 

outflow terminal speed, assumed to be 1500 km s _1 in all calculations) has been normalized 
to the average mechanical luminosity for instantaneous starbursts with a Salpeter initial 
mass function, sources between 1 and 100 Mp and with ages less than 10 Myr, Lnsb = 
3 x lO 4O (MAr5B/lO 6 M ) erg s -1 (ILeitherer et al.lll999l ). The reinserted gas was assumed to 
have a negligible angular momentum and thus the flow could be solved in spherical symmetry. 

The mass and energy deposition rates per unit volume are q m = 3Mnsb/^R%sb anc ^ 
q e = 3L^sb/4^Rnsbi respectively, where M NSB is the total mass deposition rate within the 
starburst volume (see Table 1). These values are added, at each time step, to the computed 
density and total energy in every cell, p Q u and e ta t,oid, respectively, when located inside the 
starburst volume using the following procedure: p new = p id + q m dt, and the velocity is 
corrected so that the momentum is conserved 
corrected to conserve the total energy, e^ m %d = £tot,oid — 
i s inse r ted as a form of int ernal energy e^ new = ei, m id + Qedt (see, ITenorio-Tagle et al. 
( 120071 ); I Wunsch et al.l (120081 )). The velocity of the flow at each radius is updated according 
to v new = v mid — a g dt, where a g = GM r /r 2 is the gravitational acceleration at each radius 
(see section 2). The computational domain extends over the interval (Ri n , R out ), where 
< Ri n <^ Rnsb < Rout- An open boundary condition was adopted at both ends of the 
computational grid. 

Our reference models are presented in Table 1. Here column 1 is a reference to the model, 
columns 2, 3 and 4 present the radius, mass and mechanical luminosity of the considered 
starburst, respectively. The ratio of the starburst mechanical luminosity to the critical 
luminosity and the total mass deposition rate inside the starburst volume are shown in 
columns 5 and 6. The mass of the central SMBH, unless explicitly mentioned, was assumed 



new 

, v mid = v i d p i d / p new ; the internal energy is 
" /2, and the new energy 



Pneu>V mid / 
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Table 1: The input models 



Model 




Mnsb 


log (L NS b) 


LNSB/L cr it 


Mnsb 




(PC) 


(10 8 M ) 






(M yr- 1 ) 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


la 


40 


2.0 


42.778 


0.5 


8.45 


lb 


40 


4.0 


43.079 


1.0 


16.89 


lc 


40 


6.0 


43.255 


1.5 


25.34 


Id 


40 


6.8 


43.302 


1.7 


28.72 


le 


40 


8.0 


43.380 


2.0 


33.79 


2a 


10 


0.3 


41.954 


0.27 


1.27 


2b 


10 


0.5 


42.176 


0.45 


2.11 


2c 


10 


1.0 


42.477 


0.9 


4.22 


2d 


10 


1.67 


42.698 


1.5 


7.03 


2e 


10 


2.22 


42.823 


2.0 


9.38 


2f 


10 


2.75 


42.916 


2.5 


11.61 


2g 


10 


10.0 


43.477 


9.0 


42.23 


2h 


10 


20.0 


43.778 


18.0 


84.46 



The starburst input parameters for the simulations. Column 1 is a reference to the models. The radius 
(Rnsb), mass (Mnsb), logarithm of mechanical power (Lnsb, measured in erg s _1 ), the ratio of the 
starburst mechanical power to the critical mechanical luminosity (Lnsb / Lcrit) and the total starburst mass 
deposition rate (Mnsb) are presented in columns 2 to 6, respectively. All starburst models have a central 
SMBH with a mass M BH = 10 8 M Q . 
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to be Mbh = 10 8 M in all calculations. The computational domain for models la - le 
extends radially from R in = 0.1 pc to R out = 50 pc. The inner and outer radii of the 
computational domain in the case of models 2a - 2h are 0.05 and 20 pc, respectively. 1000 
grid zones were used in all calculations. The resolution convergency was tested in the case 
of the most energetic model (2h), carried out with 1000 and 3000 grid cells resolution. 
The results with both resolutions are in excellent agreement over the whole computational 
domain. 



3. Results 

The results of the simulations are summarized in Table 2. Where, column 1 is a reference 
to the model, column 2 presents the resultant stagnation radius. The total mass deposition 
rate is shown in column 3 and should be compared with the resultant accretion rate and the 
mass outflow in the starburst wind, presented in columns 4 and 5, respectively. Column 6 
shows the SMBH luminosity normalized to the Eddington limit (L Edd = 1.3 x 10 38 M BH Mq 1 
erg s _1 ). The ram pressure of the outflow at the starburst edge (P ram = pu 2 ) is presented in 
column 7. 



3.1. Comparison of the numerical and the semi-analytic solutions 

In order to test our numerical code, several simulations were carried out for starbursts 



with Lnsb < L crit to compare them with the semi-analytic model of ISilich et al.l ( 120081 ). 
Figure 1 presents, as an example, the results of the semi-analytic (solid line) and numerical 
(open circles) calculations for case 2c in Table 1. There is a good agreement between the two 
methods, as shown in panels a - c, for the stationary run of velocity, density and temperature, 
respectively. The value of the stagnation radius, marked by the dotted line in panel (a), 
is R st = 2.8 pc, only about ~ 1.5% less than the value predicted by the semi-analytic 
model. The stationary solution shows how the matter deposited by massive stars inside the 
stagnation volume ends up falling towards the center and fuels the SMBH. On the other hand, 
matter reinserted between the stagnation radius and the starburst edge is steadily accelerated 
to reach its sound velocity at the starburst edge and then it expands supersonically forming 
the starburst wind. Note that as matter falls to the center its density grows orders of 
magnitude due to convergency alone (panel b). The temperature increases also very sharply 
(panel c) due to the violent compression induced by the rapidly in-falling matter. 
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Table 2: The predicted accretion rate and the power of the wind 



Model 


R s t 


M NS B 


M acc 




M w 


L a cc/ LEdd 


p 

± ram 




(PC) 




(M© yr- 1 ) 








(10^ 7 dyn cm- 2 ) 


(i) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


la 


2.1 


8.45 


1.88 x 10" 


-3 


8.44 


8.2 x 10" 4 


1.96 


lb 


4.7 


16.89 


4.33 x 10- 


-2 


16.86 


1.9 x 10~ 2 


3.82 


lc 


11.4 


25.34 


0.93 




24.72 


0.4 


5.39 


Id 


14.4 


28.72 


1.38 




27.33 


0.6 


5.89 


le* 


17.5 


33.79 


2.88 




30.89 


1.2 


6.57 


2a 


1.3 


1.27 


5.09 x 10" 


-3 


1.26 


2.2 x 10~ 3 


4.66 


2b 


1.7 


2.11 


1.02 x 10" 


-2 


2.07 


4.5 x 10~ 3 


7.67 


2c 


2.8 


4.22 


0.14 




4.09 


6.1 x 10~ 2 


14.64 


2d 


4.4 


7.03 


0.59 




6.45 


2.6 x lO" 1 


22.07 


2e 


5.3 


9.38 


1.36 




8.00 


0.6 


26.83 


2f* 


5.9 


11.61 


2.38 




9.23 


1.04 


30.53 


2g* 


8.2 


42.23 


23.28 




18.74 


10.2 


56.16 


2h* 


8.9 


84.46 


59.54 




25.48 


26.0 


77.43 



The predicted accretion rate and the power of the wind. Column 1 is a reference to the models. The results 
of the calculations: the value of the stagnation radius (R s t), total starburst mass deposition rate (Mnsb), 
the calculated mass accretion rate (M acc ), the rate at which matter flows away from the starburst as a super 
wind (M w ) 1 the stationary SMBH accretion luminosity normalized to the Eddington limit (L acc / LEdd) and 
the ram pressure of the wind {P ram ) are shown in columns 2 to 7, respectively. Asterisks mark those models 
for which the calculated accretion luminosity exceeds the Eddington limit. 
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Fig. 1. — Test calculations. The numerical hydrodynamic solution (circles) for model 2c 
(see Table 2) is compared with the semi-analytic results (solid lines). Panels a-c show the 
run of the stationary velocity, density and temperature across the radial direction. In panel 
(a), the dotted and dashed vertical lines mark the location of the stagnation radius and the 
nuclear starburts radius, respectively. 
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3.2. The numerical solution above the threshold line 




Fig. 2. — The comparison of the pressure gradient to the gravity force inside the starburst 
volume. Panel (a) shows the ratio of the pressure gradient to the gravity force for starbursts 
of the same mass and radius (model 2c) but with different gas metallicities. Here solid, 
dashed and dotted lines correspond to Z = Z Q , Z = 5Z Q and Z = 10Z Q , respectively. The 
intersection of the curves with the thin horizontal lines marks the position of the stagnation 
radius. The stagnation radius moves to a larger distance from the center as the cooling 
rate becomes larger. Panel (b) shows the same ratio for starburst with different energy 
deposition rates (or different masses) but the same (solar) metallicity, the solid and dashed 
lines correspond to models 2c and 2e in Table 2, respectively. 

The hydrodynamic solution for the matter reinserted by massive stars within an evolving 
young massive starburst in presence of a central SMBH is always bimodal, whether one 
considers starbursts below or above the threshold line. The main difference is that above 
the threshold line strong radiative cooling becomes the physical agent that defines where 
the stagnation radius lies. The stationary location of the stagnation radius is defined by 
the balance between the gravitational force (F g ) and the outward thermal pressure gradient 
(dP/dr), which naturally, is strongly affected by energy losses. Figure 2 presents the ratio 



- 10 - 



of the pressure gradient to the gravity force as a function of distance from the center of the 
starburst. At the stagnation radius dP/dr = F g and thus the intersection of lines which 
display this ratio with the thin horizontal line marks the position of the stagnation radius 
for various cases. Figure 2a shows the ratio of the pressure gradient to the gravity force for 
starbursts with identical mass and radii (equal to those of model 2c in Table 1) when the 
thermalized gas was assumed to have different metallicities. In these cases the displacement 
of the stagnation radius to larger and larger values is promoted by the increasingly larger 
amount of energy lost through radiative cooling within the considered starburst. Similarly, 
radiative cooling is enhanced as one considers more massive starbursts. These reinsert more 
material per unit time and thus lead to a more significant radiative cooling, as shown in 
Figure 2b, for cases 2c and 2e. In the region r < R s t, the outward pressure gradient is not 
able to compensate the gravity force and then all matter reinserted within the stagnation 
volume falls towards the central SMBH. On the other hand, in the region R st < r < Rnsb the 
pressure gradient exceeds the gravity force and hence the matter deposited there accelerates 
outwards and conf o rms a supersonic wind. Similar trends were noticed in the 3D results from 



Schartmann et al.l (120091 ) when considering the mass and energy input rate from planetary 
nebulae and type I SN in evolved clusters surrounding a SMBH. Here however, we conclude 
that the amount of matter which fuels the central SMBH and that which forms the starburst 
wind, both depend directly on the location of the stagnation radius. 

Figure 3 shows the results of numerical simulations for starbursts with a mechanical 
luminosity larger than L crit (models 2e, 2g and 2h). Here the upper, middle and lower pan- 
els present the quasi- stationary distribution of the flow velocity, density and temperature, 
respectively. As one considers more energetic (or more massive) starbursts, the larger den- 
sities (see middle panels) promote a faster radiative cooling within the thermalized plasma 
and this results in a smaller pressure gradient and thus in a further displacement of the stag- 
nation radius towards larger distances from the starburst center. This is shown by vertical 
dotted lines in the upper panels. 

The structure of the accretion flow for starbursts above the critical line presents some 
distinct features. In particular, the temperature distribution is different from that in the case 
of starburts below the threshold line. It drops smoothly within the starburst region until 
a thermal instability develops within the accretion flow. The temperature then suddenly 
drops to the minimum permitted value as shown in the bottom panels in Figure 3. 

Faster cooling leads also to a smaller wind speed. If one measures at a distance r = 
2Rnsb it is 1136 km s" 1 in the less energetic considered model (2e) and 629 km s^ 1 in the 
most energetic case (model 2h), instead of the 1500 km s" 1 expected in the adiabatic case. 
Note that for the most energetic case (model 2h), the temperature drops also suddenly to 
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Fig. 3. — The hydro dynamic solution for starbursts with a mechanical luminosity Lnsb > L cr n 
and a central SMBH. Panels from left to the right correspond to models 2e, 2g and 2h, respec- 
tively. Upper, middle and lower panels display the stationary velocity, density and temperature 
distributions, respectively. Dotted and dashed lines in the upper panels mark the location of the 
stagnation and the starburst radii, respectively. The stagnation radius appears further out for more 
energetic (and thus more massive) starbursts, because the strong radiative cooling depletes rapidly 
the temperature to the minimum allowed in the calculations, T m j„ = 10 4 K, and this results in a 
sudden lost of pressure. The radius of the thermally unstable zone also grows as the mechanical 
power of the considered starburst increases. 
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10 4 K in the free-wind region. 



Note that despite the continuous input of energy, the gas temperature remains then at 
the same value all the way, as it falls towards the center. Note also that the thermal instability 
appears at larger distances from the center as one considers a larger mass deposition rate 
or a more massive and luminous starbursts (Figure 3, bottom panels). However, in all 
cases the cold and supersonically in-falling flow is well restricted to the central regions of 
the starburst, well within the stagnation volume. This is most probably the reason why 
the semi-analytic method is able to find with great accuracy the location of the stagnation 
radius, even for star bursts above the thresho ld line. Figure 4 presents the Jeans radius, 
Rj = 0.5c s (7r/Gp) 1//2 (IClarke fc Carswelll 120071 ). where c s is the local sound speed calculated 
at each r inside the stagnation volume for our most massive models 2g and 2h. In both 
cases Rj would be larger than r if the central SMBH is able to photoionize the accretion 
flow and the gas temperature cannot drop below 10 4 K. However Rj may be smaller than 
r, and thus t he accretion flow m ay become gravitationally unstable if the temperature falls 
below 10 3 K flWada et al.l 1200a ). Note however, that for a given luminosity, the location of 
the stagnation radius is independent of the minimum temperature allowed in the flow and 
therefore, M acc does not depend on the value of T min , unless one quantifies star formation in 
the accretion flow, which is beyond the scope of the present study. Thus, the accretion flow 
is always gravitationally stable in the quasi-adiabatic regime, below the threshold line, and 
becomes progressively more unstable above the critical line as one considers starbursts with 
a larger mechanical luminosity. 

The re-inserted matter thermalized inside the starburst region may contribute to the 
observed X-ray emission. The 0.2 - 8.0 keV X-ray luminosity from the whole computational 
domain is (e.g. ISilich et al.l (120051 )): 



47T 



r 2 n e niA x {T, Z)dr 



(4) 



A' i? 



where the electron and ion number densities are n,- 



■ ,,„,., Hirm. is the average mass 

per ion, A X {Z,T) is the X-ray emissivity (see IStrickland k, Stevens! (120001 )). R m i n and R ma x 



are the inner and outer boundaries of the computational domain. This is equal to about 
2.5 xlO 41 erg s -1 and 2 x 10 42 erg s _1 for cases 2c and 2e, respectively. The luminosity of 
the infalling matter is even smaller. It is about 3.9 x 10 40 erg s -1 and 1 x 10 42 erg s _1 for 
cases 2c and 2e, respectively. This emission is orders of magnitude smaller than the SMBH 
luminosity, L acc = r] acc M acc c 2 , where r] acc = 0.1 is the accretion efficiency and c is the speed 
of light, which is: L acc « 8 x 10 44 erg s -1 in case 2c and L acc « 8 x 10 45 erg s" 1 in case 
2e. And thus, it would be hard to quantify the i nfalling matter contrib ution to the total 
X-ray emission. In this case, low luminosity AGNs (jTerashima et al.ll2002l ) seem to be better 
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Fig. 4. — The ratio of the Jeans radius to the radius of the flow as a function of distance 
to the nuclear starburst center. Panels a and b present the results of the calculations for 
models 2g and 2h, respectively. The flow is gravitationally unstable if Rj < r. The solid, 
dashed and dotted lines present \og(Rj/r) for different minimum temperatures allowed in 
the calculations: T min = 10 4 , 10 3 and 10 2 K, respectively. 
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candidates to show the infalling matter X-ray emission, as in the case of the Seyfert 2/LINER 
galaxy NGC 4303. This shows the Raymond-Smith soft X-ray emission (kT w 0.65 keV) 
originating in the core of the galaxy with r < 15 pc, coi ncident with a young ( age a round 4 
Myr, and a 3 pc radius), nuclear super star cluster (see iJimenez-Bailon et all ( 120031 )). 



4. The starburst wind and the SMBH accretion rate and luminosity 

The hydrodynamic solution discussed in the previous section allows one to calculate the 
accretion rate and thus the SMBH luminosity for each model presented in Table 1. The 
time-dependent accretion rates for models 2c - 2f calculated as the mass flux through the 
inner grid boundary are shown as examples in Figure 5. At t = Myr the accretion rate, 
M acc , is equal to zero because a stationary wind solution with R st = was used as the initial 
condition (see section 2.1). However, the accretion rate grows rapidly and reaches an average 
value of 0.14 M yr" 1 , 0.59 M yr" 1 , 1.36 M yr" 1 and 2.38 M yr" 1 for models 2c, 2d, 2e 
and 2f, respectively. Note that, the accretion rate grows due to the larger stagnation volume 
and the larger mass deposition rate from more massive starbursts. After ~ 0.1 Myr the 
total mass is conserved Mnsb = M acc + M w , see Table 2, columns 3, 4 and 5, respectively. 
Consequently, the fraction of the deposited mass expelled as superwinds from the starburst 
region decreases for more energetic starbursts although in absolute values it grows with the 
mass of the considered starburst (see Table 2). 

The stationary accretion rates onto the SMBHs, M acc , and the stationary rate at which 
matter is ejected as a starburst wind, M w , obtained through the numerical integration of 
the flow equations and normalized to the total mass deposition rate, M^sb, are shown in 
Figure 6a as a function of the normalized starburst mechanical luminosity, L^sb/ L cr i t . The 
circles represent the results from the numerical simulations of models 2a - 2e. The solid 
and dashed lines present the semi-analytic M acc and M w , respectively. The mass accretion 
rate onto the SMBH grows more rapidly when the starburst mechanical power exceeds the 
critical value, (L^sb > L crit ). This leads to a rapid increase in the central SMBH luminosity, 
L aC c = VaccM acc c 2 , as shown in Figure 6b. There, the open circles result from our numerical 
simulations, the solid and dotted lines present the semi-analytic results for starburst with 
Rnsb = 10 pc (models 2a - 2e; see Table 1) and Rnsb = 40 pc (models la - Id), respectively. 
The cross symbols mark the critical luminosity value {L^sb — L CTit ). Note, that the accretion 
rate and the SMBH luminosity obtained numerically are in a good agreement with those 
predicted by the semi-analytic model, even for starbursts with Lnsb > L crit . This implies 
that the semi-analytic calculations lead to the correct value of the stagnation radius and thus 
may be used to estimate both the starburst wind power and the accretion onto the central 
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Fig. 5. — The time evolution of the mass accretion rate. The accretion rate onto the central 
10 8 M black hole for models 2c, 2d, 2e and 2f; dotted, dash-dotted, dashed and solid lines, 
respectively. The labels indicate the L^sb/ L cr i t values. Note that after ~ 0.1 Myr, M acc 
remains almost constant and these values are reported in Table 2 together with M NSB and 
the resultant M w . Note also that, the accretion rate associated to model 2f exceeds the 
Eddington limit. The small amplitude oscillations observed, mainly in models 2e and 2f, are 
numeric artifacts. 

SMBH and its corresponding luminosity, in all cases (above and below the threshold line). 

Note that in both sets of calculations (Rnsc equal to 10 pc and 40 pc, with the assumed 
Va,oo = 1500 km s _1 ) the accretion rate reaches values ~ 1.4 M yr _1 when L NS b ~ 2L crii . 
This could result in ~50% increase in the mass of the SMBH after ~ 50 Myr. Note also 
that the calculated accretion luminosity exceeds the Eddington limit (see Table 2, models 
le, 2f, 2g, and 2h) when the starburst mechanical luminosity is just about twice its critical 
luminosity. Certainly, the accretion rate and hence the SMBH luminosity could be reduced 
if additional physics are included in the model. For example, one could think on a 2D 
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or 3D geometry that could account for the radiative and/or mechanical feedback from the 
ce ntral AGN and the redist ribution of the net angular momentum in the accretion flow (e. 
g. ISchartmann et al.l ( 120091 )). However, our ID model accounts for a realistic deposition of 
mass and energy around a central SMBH, and hence the results here presented give a good 
estimate of the accretion rate upper limit. And more important of all, the model establishes 
a direct interplay between nuclear starbursts and their central SMBHs. A direct interplay 
in which all the re-inserted matter unable to join the superwind, becomes available to the 
SMBH. The starburst wind, on the other hand, is sufficiently powerful in all to 
significantly re-structure the host galaxy ISM, leading perhaps to a thick ring, along the 
plane of the galaxy, and to a supergalactic w ind along the host galaxy symmetry axis (as in 



Tenorio-Tagle fc Munoz-Tunor] (119971 . Il998h ) 



A simple estimate of the wind power can be obtained from its ram pressure (-Pram) &t the 
starburst edge, see column 7 in Table 2. This is, in all cases, many orders of magnitude larger 
than the typical ISM pressure in our Galaxy (~ 10 ~ 12 dyn cm -2 ). It also exceeds by almost 
three orders of magnitude the pressure exerted by a one particle per cubic centimeter ISM, 
freely falling onto the starburst (Pism = Pi SM v2 ff) with Vff = [2G(M BH + M NSB ) / 'Rnsb] 1 ^ 2 '■ 
The implication is thus that the resultant winds are to lead to the build up of superbubbles 
and probably to supergalactic winds, preventing, in most cases, the falling of the ISM onto the 
nuclear starburst. Perhaps only in the case of an extremely dense ISM (pism ~ 10 -20 — 10~ 21 
g cm -3 ) freely falling onto the central starburst would modify the structure of the outflow. 



5. Conclusions 

By means of ID numerical simulations and semi-analytic estimates, we have worked 
out the stationary hydrodynamic solution for the matter reinserted by stellar winds and 
type II supernovae from a young, massive and compact starburst in presence of a central 
SMBH. The solution is bimodal in all cases, below and above the threshold line (L crit ), with 
a stagnation radius (R s t) which defines the outer boundary of the accretion flow onto the 
SMBH as well as the inner boundary of the starburst wind. 

We have shown that at the stagnation radius, the force of gravity perfectly balances 
the outward pressure gradient acquired by the thermalized reinserted matter. We have 
also shown that radiative cooling becomes an important issue for massive starbursts with a 
mechanical luminosity above the threshold line (Lnsb > L CT u). In all these cases, radiative 
cooling depletes the pressure established through thermalization of the injected matter and 
this leads to the development of a thermal instability in the accretion flow. The stagnation 
radius then moves rapidly, towards the starburst boundary, with the mass of the considered 



-17- 




-0.6 -0.4 -0.2 -0.0 0.2 0.4 41.5 42.0 42.5 43.0 43.5 

log (L N sB/L cr . it .) log (L NS b [erg s" 1 ]) 



Fig. 6. — The semi-analytic and numerical predictions for the accretion rate and the SMBH 
luminosity. Panel (a) presents the calculated mass accretion rate and the rate at which 
mass is expelled as a superwind (circles) for starbursts below and above the threshold line 
(models 2a - 2e). These are compared with the semi- analytic predictions (solid and dashed 
lines, respectively). All rates have been normalized to the total starburst mass deposition 
rate, Mnsb- Panel (b) shows the SMBH luminosity normalized to the Eddington limit. 
The circles represent results from the numerical simulations, solid and dotted lines show 
semi-analytic calculations for starburst with 10 pc (models 2a - 2e) and 40 pc (models la - 
Id), respectively. An accretion efficiency of r] acc = 0.1 was used in the calculations. Cross 
symbols represent starbursts with the critical energy (L crit ) input rate. 
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starburst. In all simulations with Lnsb > L critj strong radiative cooling occurs at a radius 
interior to the stagnation radius. Radiative cooling re-structures the inner accretion flow 
lowering the temperature to the lowest allowed value, T = T min . From then onwards and 
despite the continuous input of energy, the rapid velocity increase leading to a rapid density 
enhancement keeps and sustains the in-falling gas temperature at T = T min . It is the 
larger mass deposition rate, provided by more massive starbursts, what triggers the onset of 
strong radiative cooling and with it the shift of the stagnation radius towards the starburst 
boundary. This results in a rapid increase of the central SMBH luminosity for starbursts 
further above the critical threshold in the Lnsb — Rnsb — Mbh parameter space. 

The larger mass deposition rates provided by more massive starbursts also leads to more 
powerful starburst winds and an estimate of their mechanical power rules out the possibility 
of the ISM feeding the SMBH, at least during the type II supernova era. 

Clearly, spherically symmetric calculations, as the ones presented here, cannot account 
for the redistribution of the net angular momentum in the accretion flow. Nevertheless, 
they provide a good estimate of the upper limit to the accretion rate onto the central black 
hole, while pointing to a direct physical link between nuclear starbursts and the central 
SMBH luminosity. Our calculations do realistically account for the symmetric deposition 
of mass and energy from massive stars around the central object. This suggests that in a 
more realistic 2D or 3D geometry, able to account for the redistribution of the net angular 
momentum, the hydrodynamics would still lead to a bimodal solution with an accretion flow 
and an outward wind. In such a case however, the residual angular momentum could favor 
the formation of a gaseous disk well contained within the nuclear starburst region. We shall 
consider some of these issues in a future communication. 
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